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Abstract. We derive on the basis of Bayes' theorem a simple but general 
expression for the retrodicted premeasurement state associated with the 
result of any measurement. The retrodictive density operator is the 
normalised probability operator measure element associated with the result. 
We examine applications to quantum optical crj^tography and to the optical 
beam splitter. 
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1. Introduction 

Prediction is the determination of future events on the basis of present 
knowledge. If our knowledge is insufficient to allow a precise determination of future 
events, then we must be satisfied with predicting probabilities. Often we attempt to 
determine past events on the basis of present knowledge. A simple example is the 
reconstruction of a message from a received signal. Again, if our knowledge is 
insufficient to allow a precise determination of past events, we must be satisfied with 
retrodicting probabilities. In our communications example, this would arise if the signal 
were corrupted by a noisy channel [1]. 

In quantum mechanics there is always an element of probability. In the 
predictive formulation, even if we know the state of a system precisely before 
measurement and have full knowledge of the measurement device, the outcome of 
measurement cannot, in general, be determined with certainty. It is customary, following 
von Neumann [2], to think of the evolution of the system in terms of a deterministic, 
time-reversible, unitary part and a sudden probabilistic change upon measurement. In 
the predictive formulation we assign the premeasurement state of the system by 
prediction on the basis of the outcome of some previous preparation event. In the 
retrodictive formalism [3-6] the state of the system between preparation and 
measurement is assigned retrodictively on the basis of the outcome of the measurement. 
The sudden change occurs at the time of preparation, so the retrodictive formalism is not 
simply the time reversal of the predictive formalism. Careful use of the retrodictive 
formalism will yield the same observable correlations as those from the predictive 
formalism, so the fact that the two formalisms ascribe different states to the system 
between preparation and measurement merely emphasises that the state is no more than 
a mathematical convenience embodying knowledge about the system rather than having 
a physical existence in its own right. 

In the predictive formalism, the most general description of a measurement 
involves a probability operator measure (POM) [7]. In the retrodictive formalism, we 
can similarly describe the preparation in terms of a preparation POM. In this paper we 



use Bayes' theorem [8] to show how the retrodictive premeasurement density operator is 
simply related to the associated measurement POM element and how the (more usual) 
predictive density operator is similarly related to the preparation POM. We concentrate 
on situations in which the system does not evolve between preparation and measurement 
use the relation derived in some quantum optical applications. 

2. Bayes' theorem 

It is well known that probabilities are not objective but depend on the amount of 
knowledge available. Consider, for example, picking the winner in a ten-horse race. 
With no further knowledge, the probability of any particular horse winning is 1/10. 
Given the extra knowledge that one horse has been scratched, the probability of any one 
of the remaining horses winning is 1/9. An unscrupulous person who has fixed the race 
might be able to predict the winner with (near) unit probability based on this privileged 
additional information. Bayes' theorem [8] determines how additional knowledge 
changes the probabilities from an a priori distribution to an a posteriori distribution. It 
is important not to attach any temporal significance to the terms a priori and a 
posteriori: they simply mean with less knowledge and with more knowledge 
respectively. Prediction and retrodiction, on the other hand, do have temporal meanings. 
Prediction is determination of probabilities of later events given probabilities of earlier 
events. Retrodiction is determination of probabilities of earlier events given probabilities 
of later events. Bayes' theorem gives a means of predicting and retrodicting and also of 
connecting the two. 

Consider two sets of possible events {a,} and {bj}. Bayes' theorem states that 
the probability of two particular events a^ and bj both occurring is 



P(a,,Z7.) = P(a,\bj)P(bj) = P(bj\a,)P(a,) , 



(1) 



where P(x\y) means the conditional probability of x occurring given that y occurs. Let 
the event a precede the event b. Then the predictive and retrodictive conditional 
probabilities are given respectively by 



P{b\a^)= ^ " ^\ (2) 
' ' P(a,) 



P(a.\b.) = —^^ . (3) 
' ' P(b ) 



We usually have access to the predictive conditional probability, but can use B ayes' 
theorem to obtain the retrodictive conditional probability as 



P(a^b^ = ' ' ' = ^ ^ ' ' . (4) 

P(bj) P(bj\a,)P(a,) 

k 



In the absence of any knowledge of preference among the events a, , we assign them 
equal a priori probabilities and in this case the retrodictive conditional probability is 

P(b:\a.) 

P{a\b.)= ^ ' ■ (5) 



As an example, consider the following problem. An employee makes a choice 
each day between going to work by bus or by train. The a priori probability of the bus 
or train being chosen is P{B) or P(T) = l- P(B) respectively. The predictive 
conditional probabilities of being late if the bus or train is chosen are P{late\B) and 
P(late\T) respectively. Before the employee arrives, the boss can use these 
probabilities to predict the probability that the employee will be late as 



PQate) = P(late\B)P(B) + P(late\ T)P(T) 



(6) 



When the employee indeed arrives late, the boss can then use Bayes' theorem to 
retrodict the probability that the employee has taken the bus as 



, P(late\B)P(B) 
P(B\ late) = -^—r^^^ , (7) 
PQate) 



where the denominator can be calculated from (6). If the boss has no knowledge of 
preference in the choice of transport, then he assigns a probability of 1/2 to each. This 
leads him to calculate the probability of his employee being late to be 

Pilate) = \ P(late\ B) + \ P(late\ T) . (8) 

Then if the employee is late, the boss would retrodict the probability that the employee 
took the bus to be 



Ti/ni, ^ P(late\B) 

P(B\late) = . (9) 

P(late\B) + P(late\T) 



Naturally, the boss can change this probability by acquiring additional information for 
example by asking the employee which route was taken. 

3. Application to quantum retrodiction 

The usual, that is predictive, description of quantum mechanics involves a state 
and observables. The state is assigned on the basis of the preparation procedure. The 
state of the system then evovles until it is measured. In this paper, in order to highlight 
the roles of state preparation and measurement, we shall consider only the case where 
the system does not change between preparation and measurement. Two simple 
examples of this case are when measurement immediately follows preparation or when 
the Hamiltonian is zero. The observables are quantities associated with a measuring 
device. The probabilities associated with the results of observations are depend on 



measurement utilises a probability operator measure (POM) [7]. A measurement POM 
is characterised by its elements Ylj . These are operators with only positive or zero 

eigenvalues and have the property 

2^n, = i, (10) 

where i is the unit operator in the state space of the system. The probability that a 
measurement will result in the outcome bj given that the event a. that the system was 

prepared in the state with density operator is 

p(bj\a,) = Tr(pr%) . (11) 

For a simple von Neumann measurement, the POM elements become projectors onto 
the eigenstates of the operator corresponding to the measured observable so that this 
probability becomes the expectation value of the density operator in the associated 
eigenstate. 

Our aim is to develop a retrodictive picture corresponding to the above. Let 
Pia^) be the probability that the event occurred (corresponding in the predictive 
picture to the preparation of the state pf^). We say that the preparation device is 

unbiased if 

j;P(a,)pr'' = ^i, (12) 

where D is the dimension of the state spaced We use the concept of an unbiased source 
to simplify our derivation of the retrodictive density operator and to illustrate more 
directly the connection with the predictive formalism. Use of an unbiased source 
corresponds to no a priori information about the source. This ensures that the 



If we wish to work with an infinite dimensional space then an appropriate Hmiting procedure must be 



retrodictive state inferred on the basis of a future measurement can be based only on the 
result of that measurement. It corresponds to the situation in the predictive picture in 
which we assign states purely on the basis of earlier preparation events. Equation (12) 
allows us to define a preparation POM with elements 



= DPiaypf^ (13) 
such that 

LS, = 1. (14) 

i 

The question now is whether we can derive a retrodictive density operator p]*"' such that 
the probability that the event a^ occurred, given that the later measurement yields the 
result bj , is 

p(a,\bj) = Tr(pf%), (15) 

for all possible preparation POM elements E. . 

We can use Bayes' theorem to determine pf^ . From (1) we have 



p/ \,\ P(bjk)P{^.) ,,,, 



Substituting (11) and (15) for the predictive and retrodictive conditional probabilities we 
find 



where we have also used the definition (13). The a priori probability of the outcome bj 
is 



(18) 



from (12). Note that this is simply the expectation value of the POM element 11^. for 

-1 ^ 

the zero information state corresponding to the density operator D 1 . 



rr(prz,) 




(19) 



We require this relationship to hold for all possible E. corresponding to aU possible 
preparation POMs and hence we must set 



is the main result of our paper. An analogous relationship between the predictive 
density operator and the preparation POM element follows from our definition (13). 

We conclude this section by examining the retrodictive picture for biased 
sources, that is sources which produce states with predictive density operators Pg^^ 
prepared with a priori probabilities Pgia^) which do not satisfy equation (12). We 
introduce operators 



I retr 

i ~ 




(20) 



is the retrodictive density operator associated with the measurement outcome b. . This 



jpred 
Bi 



(21) 



analogous to the predictive POM elements, but with 



I 



(22) 



We can use Bayes' theorem and the predictive conditional probabihty (1 1) to write 



We can rewrite this in the retrodictive picture using the retrodictive state (20) as 




The lack of symmetry between this retrodictive conditional probability and the predictive 
conditional probability (11) does not reflect any intrinsic time asymmetry arising from 
quantum mechanics or from Bayes' theorem. Rather, it reflects the fact that it is 
possible, and even usual, to use biased sources of states in experiments. It is also 
possible, however, to make predictions which are biased on the basis of future 
measurements. In such cases, we would restrict our attention to only a subset of the 
possible experimental results. The predictive conditional probability would then be 



P{bi = ^F^-S ^ ' (25) 



where the sum includes only those measurement POM elements in the subset of 
interest. This situation is not uncommon in quantum optics. Consider, for example, the 
much studied phenomenon of two-photon interference in which two photons from a 
parametric down-conversion source exhibit a non-classical effect [9]. The experiment is 
usually analysed with the aid of a state involving only two-photons. The theory of the 



parametric down-converter, however, predicts a state which is a superposition of this 
state and the vacuum plus other higher photon-number states. In considering only the 
two-photon component of the state we are restricting our description to those situations 
in which the future measurements record some photocounts. 

4. Quantum optical examples 

The retrodictive formalism can be used to analyse situations involving 
preparation and measurement for a quantum system. It will provide a more natural or 
simpler description of events in some situations than in others. The retrodictive 
interpretation of phenomena may be quite different from the more familiar predictive 
interpretation but is no less valid. We have already described the retrodictive 
interpretation of photon antibunching and of the Kocher-Commins experiment [4]. 
Here, we derive the retrodictive states arising in quantum crj^tography and optical fields 
superposed by a beam splitter. 

4.1. Quantum cryptography 

A very important and obvious application of retrodiction is communication. The 
recipient of the signal (Bob) has the task of constructing the original message sent by 
the transmitter (Ahce). In order to do this. Bob uses the knowledge gained from his 
measurement together with the characteristics of the communication channel. The most 
developed technique in quantum communications is quantum cryptography or quantum 
key distribution [10], in which Alice and Bob attempt to construct a secret shared key 
which has not been leaked to any eavesdropper. 

The first protocol for quantum key distribution was devised by Bennett and 
Brassard [11]. In the Bennett-Brassard protocol Alice prepares single photons in either 
in one of two states of linear polarisation (vertical \V) or horizontal \H)) or in one of 
two states of circular polarisation (left |L) or right \R)). These states do not change on 
propagation to Bob. The states of circular polarisation can be written as superpositions 
of the states of linear polarisation and vice versa 



\L) = 2-yi\V) + i\ H)) , \R) = 2-yi\V)- i\ H)) 

\V) = 2-'/i\L) + \R)), \H) = 2-yH{\R)-\L)) . (26) 

Alice sets her preparation device so that in any given time slot of these states is 
prepared. In order to ensure security, she selects the states with equal probabilities 
thereby providing an unbiased source. Bob chooses to measure either circular or linear 
polarisation for each time slot. He chooses randomly between these two incompatible 
observables. Let us suppose that Alice chooses L for one particular photon. On the 
basis of this choice, Alice assigns to the photon the predictive state | . She can 
then use this to predict the probabilities for the outcomes of Bob's measurements 
conditioned on her knowledge of the preparation. Alice's predictive conditional 
probabilities that Bob will find each of the four possible outcomes are 

P(b^\aj) = /'(circular measured) x |(L|L)''"''|^ = i x 1 = | 
PihW) = /'(circular measured) x \{R\LY'"'\^ = i x = 
P{^by\ai^ = P( linear measured) x |(y|L)^'^^''| =jXj = \ 
PlhK) = /'(linear measured) x \{H\Lf'''f = i x ^ = ^ . (27) 

Let us suppose that Bob measures linear polarisation in a given time slot and finds 
vertical polarisation corresponding to the POM element | V){V \ . From (20) he assigns 
to the photon the corresponding retrodictive density operator Py''=|y)(y| or 
equivalently the state vector ly)'"" . He can then use this to retrodict the probabilities for 
Alice's preparation events conditioned on his knowledge of the measurement. Bob's 



retrodictive conditional probabilities that Alice prepared each of the four possible states 
are 

P(a^\by) = P(circular prepared) x V)""|^ = i x ^ = ^ 
P(aj,\by) = ^(circular prepared) x |(i?|y)""|^ = i x i = i 
P(ay\by) = P(linear prepared) x V)"""!' = i x 1 = ^ 
P(aj,\by) = /'(linear prepared) x |(//| V)""|' = i- x = . (28) 

It is straightforward to verify that the two sets of conditional probabilities satisfy Bayes' 
theorem (1) with P{a.) = j = P{bj) ■ The protocol involves a pubUc discussion between 

Alice and Bob in which some of the results are discussed. Eavesdropper activity is 
revealed by comparing the above conditional probabilities with the observed 
correlations. Alternatively eavesdropper activity is indicated if, in any time slot, Alice's 
predictive state is orthogonal to Bob's retrodictive state. 

4.2. Beam splitters 

We shall now study the retrodictive description of the optical measuring device 
depicted in Figure 1 . This consists of two photodetectors and in the output 
modes b and c of a beam splitter with a known (predictive) state in the input mode c. 
The separate POM elements for the detectors are ft^ and ft^ corresponding to 

detection events n and m. 

The POM element for the system comprising the two detectors is ft* ft^ . To 
include the beam splitter as well we change this to U Tl^ n^?7 where U is the forward 
unitary evolution operator for the beam splitter given in terms of the creation and 
annihilation operators for the modes by [12] 



U = exp 



(29) 



The reason for inserting the forward and backward unitary operators in the order we 
have done will become clear below. Finally we wish to find a POM element which 
allows us to include the field in state as part of the complete measuring device. We 
require a POM element n„„ such that when there is an incoming (predictive) state pf^ 
in mode b we obtain the probability for events n and m 



= Tr,\lJpr'kU'K^'^)-^ (30) 

where we have used the cyclic property of the trace. The reason for this requirement is 
that the second line of (30) is the probability for these events calculated in the standard 
predictive formalism using the forward evolved predictive density matrix Upf^p^ C/^ . 
For equation (30) to hold for any ""^ we require 

^nm = TrXpJJ'K^'Jj) . (31) 

From our relationship (20) the retrodictive state of the field in input mode b immediately 
before measurement, based on the measurement outcome of events n and m, is 



We shall now consider two important special cases of the above. The first 
involves modelling an inefficient photodetector by a perfectly efficient detector and a 
completely inefficient detector . The state in input mode c is the vacuum, |0)^ ^(0| . 
The individual POM elements for the detectors and are respectively \n)^ ^{n\ 



and 1^ where the latter is the unit operator for the state space of mode c. These POM 
elements are associated with recording n photocounts in and the inevitable event that 
zero counts are recorded in respectively. The retrodictive state (32) now becomes 



^retr _ ' i \ i o o \ \ c i / c (33) 

' ~Tr,(MU'\n\,{n\®iu\Ol) ' 



We find a simple expression for this retrodictive density operator by using the operator 
identity [12] 

■.oxp(-b'b):= \0)J0\ (34) 

to give 

XO\&y\M®iu\Ol = ^^MU'ib')":c^^^^ (35) 

which upon performing the unitary transformation becomes 

MU'\n),M®Km^ = ^r]"(b')"- exp(- rib'bj.b" . (36) 

Here rj = cos^ is the efficiency of imperfect detector being modelled. The 
denominator in (33) is most easily evaluated as the trace of (36) in the coherent state 
basis which gives the result 1 / 77 . Hence the retrodictive premeasurement density 
operator for mode b, conditioned on the event that the imperfect detector recorded n 
photocounts is 



^retr ^ _L;^«^lJ^tj ; exp(- T/^^^): . (37) 



This is in agreement with our previous work in which Bayes' theorem was applied to 

analyse a single imperfect photodetector [5]. 

Our second example is the measurement method of projection synthesis [13]. 

This allows us to determine the probability distribution for any physical observable 

associated with a quantised optical field mode. To keep the discussion simple, we will 

assume that both detectors are perfect so that their individual POM elements are 
and |m)^ ^("^1 • These correspond to the events n and m in which n and m 

photocounts are recorded in the detectors and respectively. We prepare the 

reference field in input mode c to be the pure state |c)^ = Co|0)^ +Ci|l)^ + L . The 
retrodictive premeasurement state (32) then becomes 

^ ^ Jc\U'^\n).,(n\®\m) (m\U\c) 

;lretr c\l I /bb\ I l/cc\ll/c /oox 
Ph — ? ~^ ^ V • \^°) 

Tr,l{c\U'\n\M®HcMU\ch) 



We can construct any retrodictive state, containing up to n + m photons, by a suitable 
choice of the reference state. Measurement by projection synthesis involves choosing 
the reference state such that the retrodictive state is an eigenstate of the observable of 
interest. As a specific example, we let the events n and m be the registering of 1 and 
photocounts respectively. From (38) we find that the retrodictive state in input mode b 
is the pure state with density operator 

_ (c; cos 0| 1), - ic; sin 0| 0) J(co cos 6> ,(1 1 + ic, sin 6> ,(0 1) 
|Co| COS a + |Cj| sm d 



Knowledge of the retrodictive state can be used for the preparation of selected 
predictive states. This can done, for example, if the field in input mode b is prepared in 
an entangled state with a field in some other mode d. Retrodicting the state in input 
mode b from the result of the measurement and projecting this onto the entangled state, 
allows us to collapse the predictive state for mode d into an unentangled field state. As 
an examole. consider the entansled state 



^> = i(|i).|o>, + /|o),|i)J 



(40) 



which can be produced by beam-splitting a single photon. Projecting the retrodictive 
pure state (39) onto the entangled predictive state (40) results in a predictive state 



, c„cos0|O) , +c, sin0|l) , 
k).= /. . " , ' (41) 



{^cj\ cos^ d + |cj|^ sin^ 



for the field in mode d. This is the retrodictive mechanism explaining the quantum 
scissors device for truncating states which we have described elsewhere [14,6]. 

5. Conclusion 

We have used Bayes' theorem to establish a relationship between the predictive 
and retrodictive formalisms of quantum mechanics. In particular we have shown that the 
premeasurement density operator for the state retrodicted on the basis of the outcome of 
a measurement is just the normalised probability operator measure (POM) element for 
the general measuring device. As well as providing a straightforward method for 
calculating the retrodictive state, it also provides a new physical interpretation of POM 
elements. Similarly the density operator for a state predicted on the basis of the 
outcome of a preparation is the normalised preparation POM associated with that 
outcome. 

We have illustrated the use of our general formula for quantum optical 
applications involving quantum cryptography and a quantum optical measuring device 
comprising two photodetectors in the output modes of a beam splitter with some known 
state in one of the two input ports. We have been able to obtain results in agreement 
with those we have previously derived by less direct techniques and have indicated how 
the retrodictive formalism can be used, by means of entanglement, in the preparation of 
selected states. 



In this paper we have restricted ourselves to the situation in which the system 
does not change between preparation and measurement. Where the evolution of the 
system density operator is unitary it is straightforward to modify our derivation of the 
retrodictive density operator to include the appropriate time evolution operator. For 
open systems, in which the evolution involves a coupling to an environment, the 
derivation is more involved. We address this problem elsewhere [15]. At a more 
fundamental level, future work might give more consideration to the role of preparation 
POMs, or preparables, which express the properties with which a system can be 
prepared, as opposed to measurement POMs, or observables, which express the 
properties which can be observed or measured. 



Acknowledgments 

We thank Howard Wiseman, Ottavia Jedrkiewicz and Rodney Loudon for 
helpful discussions and David Maguire for a careful reading of the manuscript. This 
work was supported by the U.K. Engineering and Physical Sciences Research Council, 
the Australian Research Council and the Royal Society. 



References 

[I] Shannon, C. E., and Weaver, W., 1963, The Mathematical Theory of 
Communication (Chicago: University of Illinois Press). 

[2] VON Neumann, J., 1955, Mathematical Foundations of Quantum Mechanics 

(Princeton: Princeton University Press). 
[3] Aharonov, y., Bergman, p. G., and LEBOWITZ, J. L., 1964, Phys. Rev. 134, 

B1410; Penfield, R. H., 1966, Am. J. Phys. 34, 422; AHARONOV, Y., and 

Albert, D. Z., 1984, Phys. Rev., D29, 223; AHARONOV, Y., and ALBERT, D. Z., 

1984, Phys. Rev., D29, 228; AHARONOV, Y, and VAIDMAN, L., 1991, J. Phys. 

A: Math. Gen. 24, 2315. 
[4] PeGG, D. T., and Barnett, S. M., Quantum and Semiclass. Opt. (in press). 
[5] BARNETT, S. M, PHILLIPS, L. S., and PEGG, D. T.,1998, Opt. Commun. 158, 45. 
[6] Barnett, S. M., and PeGG, D. T., Phys. Rev. A (in press). 
[7] HelSTROM, C. W., 1976, Quantum Detection and Estimation Theory (New 

York, Academic Press). 
[8] Box, G. E. P. and TlAO, G. C, 1973, Bayesian Inference in Statistical Analysis 

(Sydney: Addison-Wesley), p. 10. 
[9] HONG, C. K., OU, Z. Y., and MANDEL, L., 1987, Phys. Rev. Lett. 59, 2044. 
[10] Phoenix, S. J. D., and TOWNSEND, P. D., 1995, Contemp. Phys. 36, 165 and 

references therein. 

[II] BENNETT, C. H., and BRASSARD, G., 1984, Proc. IEEE Int. Conf on 
Computers, Systems and Signal Processing (Bangalore), pp. 175-179. 

[12] BARNETT, S. M., and RADMORE, P. M., 1997, Methods in Theoretical Quantum 

Optics (Oxford: Oxford University Press). 
[13] Barnett, S. M., and PeGG, D. T., 1996, Phys. Rev. Lett. 76, 4148; PHILLIPS, 

L. S., Barnett, S. M., and PeGG, D. T., 1998, Phys. Rev. A 58, 3259. 
[14] PEGG, D. T., PHILLIPS, L. S., and BARNETT, S. M, 1998, Phys. Rev.Lett. 81, 

1604. 



Barnett, S. M., Pegg, D. T., Jeffers, J., Jedrkiewicz, O. and Loudon, R. 
(to be published). 



Figure Caption 



Figure 1. Optical measuring device comprising a beam splitter and photodetectors in 
the two output modes b and c. A field in a known state is in input mode c. 
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